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Abstract 


The  orieniational  order  parameter  associated  with  the  second  order  nonlinear  optical 
(NLO)  coefficient  is  calculated  by  using  a  statistical  mechanics  method  that  includes  the 
angular  dependent  intermolecular  interaction  potential.  It  is  shown  that  in  addition  to  the 
alignment  of  dipoles  due  to  the  poling  field,  the  orientational  pair  correlation  (OPC)  factor 
is  also  important  in  determining  the  orientational  order  parameter.  When  OPC  is  included, 
the  order  parameter  becomes  a  function  of  the  NLO  chromorphore  concentration.  The 
concentration  dependent  part  of  the  order  parameter  is  shown  to  be  proportional  to  the 
Kirbvood  g-factor  associated  with  the  orientational  correlation  of  the  electric  dipoles.  In  the 
weak  field  limit,  it  is  shown  that  the  ratio  of  two  nonlinear  optical  coefficients,  X33/Z31. 
remains  equal  to  3,  independent  of  the  specific  nature  of  the  anisotropic  intermolecular 
potential. 
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I.  iDtroduction 


Amorphous  polymers  containing  nonlinear  optical  (NLO)  moieties  show  a  great 
promise  for  the  development  of  electro-optical  (EO)  devices  and  second  harmonic 
generation  (SHG).  However,  since  amorphous  pol>'mers  lack  the  polar  order  as  the  NLO 
moieties  in  the  polymers  are  randomly  oriented,  it  is  necessary  to  impose  polar  orientational 
order  to  break  down  the  centrosymmetry  to  induce  the  second  order  nonlinear  optical  effect. 
This  can  be  accomplished  by  poling  the  isotropic  medium  with  a  strong  electric  field.  Using 
external  field  poling,  films  that  exhibit  a  large  second  order  nonlinear  optical  effect  have 
been  prepared.^ 

In  the  usual  theoretical  description  of  the  electric  field  poling  of  the  NLO  moieties, 
it  is  commonly  assumed  that  the  dipoles  are  oriented  independently.  The  order  parameters 
associated  with  the  second  order  maaoscopic  susceptibilities  are  simply  given  by  Langevin 
functions,  L^  (a),  where  a  =  f^nEp/kT.  Here  Iq  is  the  local  field  factor  associated  with  the 
external  poling  field  Ep;  m  is  the  ground  state  dipole  moment.  One  consequence  of  the 
independent  dipole  orientation  approach  is  that  the  orientational  order  parameter, 
calculated  for  the  second  order  susceptibility,  is  practically  independent  of  the  concentration 
of  the  NLO  moieties.  Furthermore,  this  approach  necessitates  (for  the  condition  of  a<  <  1) 
the  ratio  of  the  electro-optic  coefficients,  and  that  of  the  SHG  coefficients,  d3j/d3^, 

equal  to  3,  independent  of  the  detail  of  the  poling  process.  However,  experimental  results 
obtained  for  the  EO  coefficients  r-,  and  SHG  coefficient  d-,  are  not  in  agreement  with  the 
simple  Langevin  function  description.  It  is  found  that  the  orientational  order  parameters 
depend  on  the  concentration  of  the  NLO  chromophore,^  and  the  ratio  of  the  second  order 
coefficients  may  be  greater  than  3.'^ 

In  this  paper,  we  examine  possible  causes  for  the  discrepancy.  We  use  a  generalized 
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statistical  mechanics  approach,  which,  not  only  consider  the  contribution  from  the  poling 
field,  but  also  includes  the  effect  of  orientational  dependent  intermolecular  interaction  to 
describe  the  orientational  order  parameters.  We  show  that  at  moderate  NLO  chromophore 
concentrations,  the  dipole-dipole  interaction  is  important  in  affecting  the  second  order 
nonlinear  optical  coeffidents.  Other  types  of  anisotropic  intermolecular  interactions  between 
a  NLO  moiety  and  its  sunounding  medium  are  also  considered  in  this  generalized  approach. 

In  Sec.  n  of  this  paper,  we  relate  the  second  order  molecular  hyperpolarizability  to 
the  lowest  order  maaoscopic  nonlinear  optical  susceptibility.  Two  orientational  order 
parameters,  characterizing  the  second  order  NLO  effect,  emerge  from  the  relation.  In  Sec. 
in  the  effect  of  orientational  pair  correlation  (OPC)  arising  from  the  anisotropic  (angular 
dependent)  intermolecular  interaction  is  incorporated  in  the  order  parameter  expressioa 
In  this  section,  the  specific  dipole-dipole  interaction  is  considered,  and  the  order  parameters 
are  related  to  the  Kirkwood  g-factor.  In  sec.  IV,  the  general  type  of  anisotropic 
intermolecular  interaction  potentials  is  included,  and  we  show  that  in  the  weak  poling  field 
limit  (with  a  <  <  1),  the  orientational  pair  correlation  gives  rise  to  a  pronounced  effert  on 
the  order  parameter;  however,  the  X3S/Z31  =  3  relationship  is  not  changed  by  the  presence 
of  OPC.  A  summary  and  conclusion  of  the  results  of  this  work  are  presented  in  this  section. 

II.  Second  Order  Hyperpolarizability  and  Macroscopic  Susceptibility 

We  restrict  the  discussion  to  the  second  order  nonlinear  optical  (NLO)  process.  Tb: 
origin  of  the  second  order  nonlinear  optical  effect  arises  from  the  molecular 
hyperpolarizability,  In  the  molecular  coordinate  system,  the  nonlinear  optical  response 
of  the  dipxjle  moment  of  a  representative  molecule  o,  subject  to  external  electromagnetic 
fields  at  frequencies  Oj  and  ^2,  is  given  by: 
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(tOj  +  (1)2)  ~  '^2'  *^1*  ^2)  (^1)  ^K  (''^2) 


(1) 


where  I,  J,  and  K  refer  to  the  principal  axes  of  molecule  o  with  respect  to  the  molecular 
coordinate  system.  The  h>perpolarizabLlity,  is  a  function  of  frequencies,  0)2,  and  + 
0)2-  The  frequency  dependence  of  the  hyperpolarizability  has  to  be  calculated  by  using 
quantum  mechanics.  For  SHG,  02=02“  for  the  linear  EO  effect,  we  have  o^  = 

o  and  02  =  0.  Fj  (o^)  and  Fjr  (02)  are  the  Jth  and  Kth  components  of  the  local  field  F  seen 
by  molecules  at  frequencies  o^  and  o,,  respectively.  In  the  Lorentz  model,  the  local  field  F 
(o)  is  related  to  external  field  E  (o)  by,  F  (o)  =  f^  E  (o),  where  f^  is  t‘3e  local  field  factor 
given  by  =  (n,^^  +  2)/3. 

The  corresponding  nonlinear  optical  polarization  induced  by  '.he  poling  field  and 
optical  electromagnetic  fields  is  given  by 

PCOj  +  Oj)  =  1  £  |JI<“>  (o,  +  o^  ^2) 


where  V  is  the  macroscopic  volume. 

It  is  customary  to  write  the  optical  p>olarization  in  terms  of  the  second  order 
maaoscopic  susceptibihty  as 

Pi(o,+Wj)  =  xS  (-Ui-Qj;o,,W2)  E2(Q,)Ek(Q2)  (3) 

In  Eq.  (3),  we  assume  that  the  optical  field  at  frequency  Oj  is  polarized  along  the  laboratory 
j-axis,  and  that  at  frequency  W2  is  polarized  along  the  laboratory  k-axis.  One  notes  that  in 
Eq.  (1)  the  local  field  is  used  to  define  the  molecular  byperpolarizability,  whereas  in  Eq.  (3) 


Page  5 


the  macroscopic  external  optical  fields  are  used  to  defioe  the  second  order  optical 
susceptibility.  For  SHG,  j  and  k  are  identical;  in  the  ensuing  calculation,  we  set  the  z-axis 
of  the  laboratory  coordinate  system  to  be  along  the  dc  poling  field  axis,  and  the  y-axis  to  be 
along  the  direaion  perpendicular  to  the  plane  of  incidence.  In  Eq.  (3),  the  component  i,  j 
or  k  takes  on  one  of  the  laboratory  coordinate  axes. 

The  molecular  hyperpolarizability  that  appears  in  Eq.  (1)  is  a  third  rank  tensor. 
However,  for  simplicity  we  assume,  in  this  paper,  that  only  one  major  component  of 
hyperpolarizability  along  the  3-axis  of  the  molecular  principal  coordinate  system  dominates. 
We  define  this  component  to  be  ^333  =  Under  this  simplifying  approximation,  the  second 
order  susceptibility  tensor  can  be  related  to  the  dominant  molecular 

hyperpolarizability,  /3,  in  a  simple  expression.  By  substituting  Eq.  (1)  into  (2),  and  then 
replacing  the  summation  over  all  molecules  by  an  average  over  a  statistical  distribution 
function,  we  obtain 


where  p  is  the  number  density  of  the  NLO  moieties,  aj3^^^  is  the  projection  of  a  unit  vector 
along  the  molecular  principal  axis  (the  3-axis)  of  the  representative  molecule  1  onto  the 
laboratory  axis  i.  The  angular  brackets  denote  an  average  over  the  one-particle  distribution 
function  and  is  given  by 

<a;;^  =  / 4'^  4^  4'  G({l),Ep)d(l)  (5) 

In  Eq.  (5)  and  in  the  remaining  part  of  this  pap)er,  we  shall  employ  the  symbol  {a}  to 
represent  3  Cartesian  coordinatesj^  =  (x„,  y*.  z^)  and  2  angular  coordinates  (0^,  to 
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specify  the  position  and  orientation  of  the  principal  axis  of  molecule  a.  Thus,  using  this 
notation,  we  have  Jdlaj  s  J  d4i,JJ  sin0^d6,  . 

The  one  particle  distribution  function  G({1},  Ep)  entering  into  Eq.  (5)  is  derived 
from  the  N-body  Gibbsian  distribution  function  by  integrating  over  the  rest  of  (N-1) 
molecules  with  respect  to  their  spatial  and  angular  coordinates; 

GaH.Ep)  =  ^/e  -  d{2}d{3}-d{N}  (6) 

where  21^  is  the  configuration  integral  given  by 

■  f  (7) 

Here  Uj^j  is  the  potential  energy  of  interaction  of  N  molecules.  The  potential  energy  of 
interaction  depends  on  the  positional  as  well  as  the  orientational  coordinates  of  the  N- 
molecules. 

By  poling  the  isotropic  NLO  medium  with  an  external  electric  field  Ep,  the  system 
becomes  uniaxial  and  adopts  a  symmetry.  While  in  the  general  case  there  are  18  second 
order  NLO  coefficients^  there  are  only  two  independent  nonvanishing  elements  for 
the  system,.  They  are  given  by: 
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xSx  =  xS  =  Ppfi-fi 

Equation  (8a)  is  obtained  by  setting  =  cosd^,  in  Eq.  (5)  and  Eq.  (8b) 

is  obtained  by  setting  a^^^^  =  cos  6^,  and  3^3^^^  =  a^^^^  =  sin  0j^  cos  <Pi  and  then 
performing  an  average  over  <pj^  in  Eq.  (5).  The  second  identities  introduced  in  Eqs.  (8a)  and 
(8b)  are  the  result  of  Voigt’s  contracted  notation  in  which  one  sets  the  last  two  subscripts 
zz  =  3  and  xx  =  1,  at  the  same  time  naming  the  first  subscript  z  as  3.  The  subscript  1, 
associated  with  0j^  refers  to  the  representative  KLO  chromophore  1.  Clearly,  in  the  absence 
of  the  poling  field  Ep,  the  angular  averages  over  cos  0^  and  cos^  0^  vanish  identically  due 
to  the  fact  that  the  medium  is  macroscopically  homogeneous  when  Ep  =  0. 

Generally  in  all  previous  work  on  the  order  parameter,  it  is  assumed  that  the  dipoles 
{ifi>  ■  ■  those  associated  with  the  NLO  moieties)  are  surrounded  by  a 

homogeneous,  thermodynamically  independent  isotropic  medium  so  that  the  dipoles  are  in 
thermal  equilibrium  with  the  external  poling  field  Ep.  In  this  assumption,  the  potential 
energy  of  interaction  of  the  molecules  where  the  dipoles  are  located  do  not  play  a  role,  and 
the  orientation  of  the  dipoles  is  completely  dictated  by  the  external  poling  field.  As  a  result 
,  the  intermolecular  interaction  term  in  Eq.  (6)  factors  out  of  the  N*body  distribution 
function  and  the  evaluation  of  the  order  parameters  given  by  Eq.  (5)  reduces  to  a  simple 
independent  one-body  problem.  In  this  exceedingly  simplified  case  and  X3i^^^  are 
readily  found  to  be 

a) 

— iH— .  l„=(4*l)cotha-(-i  *  ^)=L,  (9a) 
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(9b) 


^  a'  a  » 

where  Lj,  (a)  is  the  Langevin  function  of  order  n  and  is  given  by 

r’'e*“*®‘cos“0.  sin6,d6, 

L  (a)  =  ^ - ,  wiiha=fopEp/kT  (10) 

j"J'e*“*®‘  sm0jd0, 


Dearly,  for  a<  <1,  Equations  (9a)  and  (9b)  reduce  to 


(11a) 


L3:  = 


a 

15 


+  ••• 


(11b) 


Thus,  in  this  simplified  case,  the  ratio  of  to  X3i^^^  is  thus  equal  to  3.  Furthermore,  L33 
and  L3j  display  only  a  negligible  chromophore  density  dependence  through  the  local  field 
factor  f^ 


III.  EfTects  of  Orientational  Pair  Correlation 
To  obtain  an  NLO  polymer,  one  either  disperses  chromorphores  exhibiting  a  large 
hyperpolarizability  )3  in  an  isotropic  amorphous  polymer  matrix  to  form  a  guest/host  system 
or  incorporates  chromorphores  chemically  into  polymer  chains.  Since  the  NLO 
chromophores  are  also  moleailes  with  large  permanent  dipole  moments,  the  NLO  polymer 
system  is  an  assembly  of  molecules  having  large  dipole  moments.  For  this  reason,  the  dipole- 
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dipole  interaction  needs  to  be  included  in  the  analysis.  Consider,  for  example,  a  molecule 
with  a  permanent  dipole  moment  equal  to  5  debye  (a  moderate  sized  dipole  moment  for 
an  NLO  chromophore)  subjected  to  a  poling  field  of  about  lO^V/cm.  This  is  a  very  large 
poling  field  which  can  be  achieved  only  by  careful  sample  preparation  with  iom’c  impurities 
removed  the  voltage  avoid  breakdown.  For  this  static  field,  the  dipole-electric  field 
interaction  energy  is  equal  to  1.67  x  10’^^  erg.  If  the  concentration  of  the  NLO  molecules 
is  taken  to  be  p  =  10^^  molecules/cm^  and  thus  the  average  intermolecular  distance  is 
about  =  10'^  cm,  then  the  dipole-dipole  interaction  energy  would  be  about  25  x  10'^^ 
erg,  which  is  larger  than  the  electric  field-dipole  interaction  energy.  Therefore,  in  contrast 
to  the  usual  approach  by  assuming  independent  dipoles  in  calculating  the  orientational  order 
parameters,  the  intermolecular  dipole-dipole  interaction  energy  cannot  be  ignored.  In 
addition,  it  is  well  know  that  the  dipole-dipole  interaction  is  long-range.  The  long-range 
nature  of  the  interaction  is  expected  to  give  rise  to  a  significant  orientational  pair 
correlation  (OPC)  effect  in  the  second  order  nonlinear  optical  coefficient  for  an  assembly 
of  polar  molectiles.  However,  the  effect  of  the  long  range  OPC  has  yet  to  be  considered. 
In  this  section  we  first  incorporate  the  OPC  factor  into  the  expressions  of  orientational  order 
parameters  L33  and  We  consider  in  this  paper  only  the  case  in  which  the  interaction 
energy  of  dipoles  with  the  poling  field  is  small  compared  with  the  thermal  energy,  which 
allows  the  a<  <  1  approximation  to  be  used.  At  100°C,  the  thermal  energy  kT  is  5.15  x 
erg.  The  other  case  in  which  the  jx)ling  field  is  large  compared  with  kT  is  more  involved  and 
will  be  reported  in  another  publication. 

By  assuming  that  the  interaction  energy  between  the  poling  field  and  the  molecular 
dipoles  is  small  compared  to  the  thermal  energy,  we  can  expand  the  exponential  factor  in 
Eqs.  (6)  and  (7)  in  a  power  series  of  the  poling  field  keeping  only  the  terms  linear  in 
Epas 
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(H.  ipXoj 

■  E  (a,-§p)f.*-  t^c-u^kD  (12) 

Jci  l-l 

i  ^ 

By  substituting  Eq.  (12)  into  Eq.  (7),  we  find  that 

=  /e'''^d{l}-d{N}  ^  0(eJ) 

=  Z«>.0(EJ, 

where  is  the  configuration  integral  of  the  N-body  system  without  the  poling  fip’d. 
Thus,  in  the  a<  <  1  linuq  is  equal  to  the  correction  being  2nd  order  in  Ep.  As  a 
result,  Eq.  (6)  can  be  written  as 

G  ({ll.El  =  _!_(l+acos0i)  +  (N-1)  /  -  /  cos8jP2(1.2)d{2)  (14) 

^  V(47t)^  V2(4it)^  •' 

where  P2  (1,  2)  is  a  two  particle  distribution  function  which  is  obtained  by  integrating  e' 
over  the  (N-2)  molecules  with  respect  to  their  positional  and  orientational 
coordinates.  The  second  term  on  the  right  hand  side  of  Eq.  (14)  has  a  factor  of  (N-1) 
because  there  are  (N-1)  molecular  neighbors  for  the  representative  molecule  1.  The  two 
particle  distribution  is  averaged  over  the  configuration  of  the  second  particle  described  by 
{2}.  This  term  gives  rise  to  OPC  and  cannot  be  neglected. 

The  normalization  condition  for  P2  (1,  2)  is 

— ^  /  pj(l,2)dtl)  d{2}  =  1  (15) 

V^(47:)*  •' 

Thus  p2  (1,  2)  is  a  dimensionless  probability  function,  which  at  low  density  is  given  by^ 
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(16) 


pO  (12)  = 

where  u  (1,  2)  is  the  angular  dependent  interaction  potential  between  molecules  1  and  2. 

The  first  term  on  the  light  hand  side  of  Eq.  (14)  is  the  usual  independent  one- 
particle  term,  which  upon  substituting  into  Eq.  (5)  gives  immediately  Eqs.  (11a)  and  (11b) 
after  performing  a  simple  integration.  The  second  term  is  the  contribution  due  to  the  pair 
correlation  of  two  interaaing  dipoles.  We  shall  show  that  if  the  intermolecular  interaction 
i5  anisotropic,  then  this  term  can  make  an  important  contribution  to  the  second  order 
optical  coefficients. 

Substitution  of  Eq.  (14)  into  Eq.  (5)  (or  specifically  into  Eqs.  (8a)  and  (8b))  gives  the 
expression  for  the  order  parameter  by 

<cos*0,>  =  +  — ^p/  /  cos»0iCOS02  g(1.2)  d{l}d{2}  (17) 

n+2  V(47c)^ 

with  n  =  1  or  3. 

Here  g  (1,2)  is  the  orientational  dependent  radial  distribution  function,  which  is 
related  to  the  two  particle  distribution  function  P2  (1,2)  by  g(l,2)  =  p2  (1,2)  -  1.  The 
distribution  function  vanishes  when  two  particles  are  separated  at  a  large  distance.  More 
important  is  the  fact  that  in  order  to  contribute  to  the  order  parameter  <cos“  0j>,  g  (1,2) 
must  depend  on  orientational  variables.  Thus,  the  effect  of  the  intermolecular  interaction 
on  the  second  order  NLO  coefficients  lies  in  the  orientational  dependent  part  of  the  g  (1,2) 
function,  which  we  shall  call  the  molecular  pair  correlation  (MPC)  functiom  The 
orientational  pair  correlation  (OPC)  can  be  projected  out  from  the  molecular  pair 
correlation  function.  In  accordance  with  Eqs.  (8a),  (8b)  and  (17),  the  NLO  coefficients  are 
then  completely  determined  once  the  OPC  is  known.  Thus,  in  addition  to  the  poling  field 
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which  induces  orientation  of  the  electric  dipwles,  the  presence  of  OPC  also  assists  the 
orientational  ordering  of  the  NLO  chromophores.  As  shown  in  recent  experiments  in  our 
laboratory,  even  at  a  moderate  NLO  chromophore  concentration  (above  10^  molecules  per 
cm^),  the  OPC  assisted  contribution  to  the  order  parameter  cannot  be  ignored.^ 

It  should  be  pointed  out  that  the  OPC  factor  is  the  same  type  that  is  involved  in  the 
theory  of  the  static  dielectric  constant  of  polar  fluids.  The  study  of  g  (1,2)  due  to  the 
orientational  dependent  interaction  in  a  homogeneous  fluid  consisting  of  anisotropic 
molecules  has  been  formally  considered  by  several  authors.®*^®  Wartheim  has  obtained  the 
solution  of  g(l,2)  using  the  mean  spherical  qjproximation  for  the  hard-sphere  dipolar 
Uquid.^^  Several  theories  based  upon  the  cluster  expansion  by  Stell,^  and  Stell  and  Weiss^ 
have  also  been  developed  to  treat  the  hard  sphere  dipolar  fluids.  Recently,  Monte  Carlo 
(MC)  calculations  have  also  been  used  to  assess  the  accuracy  of  these  analytical  theories.^^ 
To  obtain  the  molecular  pair  correlation  (MPC)  function  one  employs  the  Omstein- 
Zemilce  equation^; 

g(U)  =  C(I,2)  4  ±  f  C(l,3)g(3,2)d(3)  (18) 

where  n  =  4r  for  a  linear  molecule,  and  is  equal  to  for  a  non-linear  molecule.  In 
Eq.(18),  C(l,2)  is  the  molecular  direct  correlation  function. 

The  solution  of  Eq.(  18)  for  a  polar  fluid  for  which  the  intermolecular  potential  is 
given  by 


(5 


u(i;i)  =  u(rij)  + 


trtj  A  ‘^12) 


‘12 


r’n 


(19) 


has  been  considereti  e)ctensively  in  the  literature^.  In  Eq.  (19),  u(rj2)  is  assumed  to  be 
spherically  symmetric,  and  the  angular  dependent  term  is  represented  by  the  dipole-dipole 
interaction. 

For  this  t>'pe  of  the  intermolecuJar  interaction  potential,  it  was  first  shown  by 
Wertheim  that  the  solution  for  g(l,2)  for  Eq.  (18)  can  be  assumed  to  be  of  the  form^^ 

g(l,2)  =  g,(ri^  +  h^  (T^  A(l,2)  +  ho  D(l,2)  (20) 

where  gs(ri2)  is  spherically  symmetric,  depending  only  on  the  amplitude  ofjj^-  ^(^>2)  = 
\  and  D(l,2)  =  3  (s^  •  r^  (S2  •  r^^)  -  (s^  •  S2).  Here  Sj  and  S2  are  unit  vectors  parallel 
to  the  dipole  moments  of  molecules  1  and  2,  respectively. 

Substituting  Eq.  (20)  in  Eq.  (16)  and  after  carrying  out  some  extensive  algebra,  one 
obtains 

<COS0J>  =  ^ 

<cos^0j>  =  ^  +  ap(~)|"h4(r)r^dr  (20b) 


One  notes  that  the  h^^'^  part  in  the  expansion  of  g(l,2)  makes  a  contribution  to  the  order 


parameters.  Thus,  even  including  the  OPC  effect,  we  obtain  in  the  weak  poling  field  limit 
rather  simple  expressions  for  and  X3i^^^  as 

Xm  = 

and 

%g  =  pPa(j-f0.)  Plb) 


where  is  the  cluster  integral  corresponding  to  given  by 

Ga  =  Y 


(22) 


We  need  not  go  into  the  discussion  of  detailed  calculations  of  h^  (r)  using  various 
theoretical  models,  but  it  suffices  to  mention  that  the  Kirkwood  g-factor  defined  by 


Sk"' 


>(M)^ 


(23) 


is  related  to  by^^ 


gjr  =  1  +  p  G^  (24) 

In  Eq.  (23),  M  is  the  total  dipole  moment  of  the  sample. 

One  may  also  include  the  dipole-dipole  interaction  between  the  NLO  chroraophores 
and  the  segments  of  the  polymer  host  in  the  calculation.  The  result  will  add  another  term 
to  Eq.  (21)  in  the  form  of  p  'G  ,  where  the  prime  designates  the  segmental  dipoles  of  the 


host  However,  the  total  OPC  contributions  from  both  the  chromophore  and  the  host  can 
be  incorporated  into  the  Kirkwood  g-factor  by  defining  Eq.(24)  as  gj^  =  1  +  p  G^+  p  'G'^. 
As  a  result  in  terms  of  the  Kirkwood  g-factor,  we  obtain  very  simple  expressions  for 
and  as 

=  pl3a  Ek/IS  (25a) 

and 

=  P^a  gj^/S  (25b) 

The  ratio  of  to  ^  ^ identical  to  the  independent  dipole  case.  Thus, 
the  inclusion  of  the  dipole-dipole  interaaion  in  the  calculation  of  the  order  parameters  does 
not  change  the  1/3  relationship.  However,  the  effect  of  OPC  gives  rise  to  a  complex 
chromophore  density  dependence  in  both  and  z^^^33.  as  manifested  in  the 

multiplication  of  the  independent  dipole  result  by  the  Kirkwood  g-factor.  Since  the 
Kirkwood  g-factor  is  related  to  the  static  dielectric  constant  c  by^ 

-  g,y  (26) 

9e 

where  y  is  a  molecular  parameter  defined  as  y  =  l5iLP  ,  Equations  (25a)  and  (25b),  in 


combination  with  Eq.  (26)  provide  a  useful  means  for  characterizing  the  molecular 
hyperpolarizability  p,  once  z^^^3i  (or  Z^^^33).  ground  state  dipole  moment  n,  and  the 
static  dielectric  constant  e  are  experimentally  determined. 

rv.  General  Anisotropic  Intermolecular  Interactions 
We  next  show  that  the  1/3  ratio  is  not  limited  to  the  dipole-dipole  interaction 
mechanism.  Rather,  it  is  the  result  of  the  weak  field  approximation.  To  show  this,  we 


Pao^ 


consider  a  general  anisotropic  intermolecular  interaction  potential.  The  anisotropic  potential 
could  either  be  short  or  long  range  or  both.  To  evaluate  the  OPC  effect  on  the  order 
parameters,  one  may  use  the  invariant  expansion  of  the  angular  dependent  radial 
distribution  funaion.^’^^  This  is  often  used  in  equilibrium  statistical  mechanics  of  complex 
fluids.  The  spherical  invariant  expansion  of  g  (1,2)  for  linear  molecules  can  be  written  as^^ 

(27) 

A 

where  r^  indicates  the  orientation  of  the  vector  -I2  ‘-Ii-  amplitude  of_r^  is  r^ 
=  hul-  The  spherical  invariant  function  is  given  by 


1.  2-12^  2li  +  r  '2«J  +  1  '2«3  +  r 

Y,...(ni)Y^^  (Q2)  Y.,^  (f,2) 


4jc 


4it 


*2  ^3) 
“2  “3 


) 


(28 


where  the  quantity  enclosed  with  circular  brackets  is  a  Wigner  3-j  symbol,  and  is 

an  orientational  independent  function. 

Substituting  Eqs.  (27)  and  (28)  into  Eq.  (16),  and  after  integrating  over  the  angular 
coordinates  we  obtain  the  expression  for  <cos  8j>  as 

<cos0,>  =  "I  +  apG,jt/3  (29) 


and 


n _ _  i  *7 


(30) 


L,,  =  W0j>  =  "I  +  pG„</5 

Hence, 

Lj,  =  ^(<cos0i>  -  W0,>J=  ^  +  I  apG,,^  (31) 

where  ^  cluster  integral  given  by 

G„o  =  ^  r  (32) 

3v/3 

which,  for  the  dipolar  liquid  is  equal  to  given  in  Eq.  (22). 

Thus,  in  the  weak  field  limit,  only  the  type  of  orientational  dependent  intermolecular 
potentials  between  molecules,  whether  it  arises  fi-om  short-  or  long-range  interactions,  that 
give  rise  to  the  dipolar  order  G^^g  contribute  to  the  second  order  nonlinear  optical 
coefficients.  The  quadruple  (i  =  2)  or  other  types  of  resultant  order  will  not  contribute  to 
the  second  order  NIX)  coefficients.  This  is  a  general  result,  independent  of  the  assumption 
of  the  nature  of  the  intermolecular  potential.  One  also  notes  from  Eqs.  (19)  and  (20)  that 
even  in  the  presence  of  OPC,  the  ratio  remains  equal  to  3.  The  order 

parameters  L33  and  L31  depend,  however,  on  the  chromorphore  density  through  the  cluster 
integral  G^^g. 

Considering  the  fact  that  <cos^  0j>  has  a  different  rotational  symmetry  from  that 
of  <cos  0j>,  one  may  wonder  about  the  invariance  result  of  the  X3/^Vx33^^^  ratio  in  the 
presence  of  the  anisotropic  intermolecular  interaction  potential.  This  situation  can  be  made 
clearer  if  one  expresses  the  order  parameters  in  terms  of  the  average  over  Legendre 
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polynomials,  Pj  (cos  Sj)  and  P3  (cos  0j).  In  terms  of  the  Legendre  polynomials,  one  can 
write  the  order  parameters  as 


L3,  =  I  {<Pj(cos0,)>  -  <P,(cos0,)>}  (33a) 

and 

L33  =  j  {3<P,(cos0i)>  +2  <P,(cos0i)>}  (33b) 


In  the  weak  field,  the  OPC  caimot  connect  the  <P3  (cos  0i)>  term,  and  as  a  result  the 
L31/L33  ratio  is  1/3,  independent  of  the  specific  nature  of  the  anisotropic  intermolecular 
interaction. 

Two  possible  mechanisms  will  bring  in  the  <P3>  term  and  results  in  the  deparmre 
from  the  1/3  ratio.  One  is  the  deviation  from  the  phase  matching  condition  that  makes  the 
second  order  susceptibility  wave  vector  dependent,  and  the  other  arises  from  a  high  poling 
field. 

The  calculations  carried  out  in  this  work  are  based  upon  a  general  assumption  of 
perfect  phase  matching  in  the  second  order  NLO  effect  As  a  result  of  this  assumption,  the 
second  order  NLO  susceptibility  does  not  dependent  on  the  wave  vector.  In  the  situation 
where  the  phase  matching  condition  is  not  met,  and  the  sample  is  spatially  inhomogeneous, 
the  second  order  susceptibility  would  become  wave  vector  dependent  In  this  situation,  the 
efficiency  of  second  harmonic  generation  will  be  low;  however,  the  X33^^Vx3/^^  =  3 
relationship  will  break  down  even  for  the  weak  field  dipolar  system.  The  other  situation  that 
will  give  rise  to  the  deviation  fi-ora  1/3  for  the  X3i^^Vx33^^^  ratio  is  the  case  of  high  poling 
field.  For  example,  when  a  =  2,  a  quick  calculation  would  show,  for  independent  dipole 
orientation,  the  X33^^Vx3i^^^  ratio  to  be  equal  to  3.5.  Intermolecular  interactions  among 
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NLO  chromophores  and  betv^ecn  ihe  NLO  chromophores  and  segmentaJ  dipoles  of  the  host 
cause  further  departure  and  result  in  an  even  larger  ratio.  We  shall  report  elsewhile  the 
calculation  of  the  high  poling  field  case. 

In  summary,  we  have  sho\^-n  that  the  orientational  order  parameters  associated  with 
second  order  nonlinear  optical  susceptibilities  depend  on  the  NXO  chromorpbore 
concentration.  This  is  due  to  the  presence  of  anisotropic  intermolecular  interactions.  The 
concentration  dependent  part  of  the  order  parameters  are  shown  to  be  related  to  the 
Kirkwood  g-factor  associated  with  the  correlation  of  electric  dipoles  in  polar  fluids.  In  the 
weak  field  limit,  the  X3-,^‘Vx3/^^  ratio  is  not  changed  by  the  presence  of  the  anisotropic 
intermolecular  interactions. 
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